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^ ! ABSTRACT 

Q , The Swift mission has discovered an intriguing feature of Gamma-Ray Burst 

\^ • (GRBs) afterglows, a phase of shallow decline of the flux in the X-ray and optical 

' lightcurves. This behaviour is typically attributed to energy injection into the burst 

d . ejecta. At some point this phase ends, resulting in a break in the lightcurve, which is 

commonly interpreted as the cessation of the energy injection. In a few cases, how- 
ever, while breaks in the X-ray lightcurve are observed, optical emission continues its 
^ . slow flux decline. This behaviour suggests a more complex scenario. In this paper, we 

' present a model that invokes a double component outflow, in which narrowly coUi- 

QQ I mated ejecta are responsible for the X-ray emission while a broad outflow is responsible 

\^ I for the optical emission. The narrow component can produce a jet break in the X-ray 

' lightcurve at relatively early times, while the optical emission does not break due to 

its lower degree of coUimation. In our model both components are subject to energy 
■ injection for the whole duration of the follow-up observations. We apply this model to 

' GRBs with chromatic breaks, and we show how it might change the interpretation of 

. the GRBs canonical lightcurve. We also study our model from a theoretical point of 

view, investigating the possible configurations of frequencies and the values of GRB 
physical parameters allowed in our model. 
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1 INTRODUCTION 

Since its launch, the Swift mission (Gehrels et al. 2004) 
has allowed us to observe the emission from Gamma-Ray 
Burst (GRB) afterglows in the X-ray and UV/Optical from 
as early as ~1 minute after the burst trigger by means 
of the X-ray Telescope (XRT, Burrows et al. 2004) and 
UV/Optical Telescope (UVOT, Roming et al. 2005). This 
unprecedented response time has allowed us to unveil the 
early behaviour of GRB afterglow lightcurves, which turn 
out to be more complex than expected. Typically, at the end 
of the prompt emission the X-ray flux F exhibits a rapid 
decay. This can be modeled with a powerlaw F ^ t~"^ with 
slope ai ~ 3 — 5 (Tagliaferri et al. 2005). This phase, which 
usually lasts hundreds of seconds, is widely interpreted as 
the tail of the prompt emission (Kumar & Panaitescu 2000; 
for a review, see Zhang et al. 2006). After that, the X-ray 
flux decays in a much shallower way, forming a "plateau" 
with a slope 02 ~ 0.1 — 0.8. The spectrum in this phase 



can be different from that observed during the fast decay, 
which indicates a different physical origin. The duration of 
the slow decline is a few thousands of seconds (O'Brien et 
al. 2006, Willingale et al. 2007). After this time, a break 
occurs and the lightcurve becomes steeper, with a powerlaw 
slope of Qf3 ~ 1.3. Indeed, this latter phase was studied well 
prior to the launch of Swift (e.g. De Pasquale et al. 2006; 
Gendre et al. 2006) and it is understood to be emission from 
synchrotron radiation, resulting from a shock produced 
by the expansion of the burst ejecta into the circumburst 
medium (Meszaros & Rees 1997). Occasionally, a further 
break may occur a few days after the trigger, leading to a 
segment with decay slope of 04 ~ 2 . This steep decay can 
be interpreted as the signature of coUimated outflow (Sari 
et al. 1999). Overall, this evolution of the X-ray flux is now 
referred as the "canonical" X-ray lightcurve (Nousek et al. 
2006). In the optical band, the flux decays with a similar 
range of slopes to those of the X-ray, with the exception of 
the initial fast decay phase, which is usually absent (Oates 
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et al. in preparation). 

The slow decay is probably the most perplexing among 
the novel aspects discovered by Swift, and several models 
have been proposed to explain it (see e.g. Zhang 2007 for 
a complete review). These models in general fall into three 
main classes: i) energy injection into the burst ejecta, either 
in the form of Poynting flux or late time shells of jecta; ii) a 
non uniform angular energy distribution in the jet or a jet 
seen off-axis, so that a fraction of the early afterglow emis- 
sion is not fully beamed towards the observer; iii) a change 
of the microphysical parameters that leads to an increase in 
the conversion efflciency of the ejecta energy to radiation. 

Puzzlingly, in a few Swift GRBs the slow decline phase 
ends with a "chromatic break" (Panaitescu et al. 2Q06a; see 
also Melandri et al. 2008): i.e. a transition from the shal- 
low to the normal decay appears in the X-ray band but is 
absent in the optical band, where the flux continues to de- 
cline at a slow rate. This feature is very hard to explain 
with any model that predicts a single origin for the X-ray 
and optical emission. In the attempt to solve this problem, 
Ghisellini et al. (2007) suggested a model in which the opti- 
cal emission is caused by the interaction between the ejecta 
and the circumburst medium, while the X-ray radiation is 
produced by internal shocks occurring in coUimated shells 
emitted by the GRB central engine at relatively late times. 
If the Lorentz factor F of these shells decreases with time, 
a "jet-like" break will be detected (in the X-ray band only) 
at the time in which F"^ = d, where 9 is the opening an- 
gle of ejecta. An alternative scenario, proposed by Genet et 
al. (2007) and Uhm & Beloborodov (2007), assumes that 
both the X-ray and optical emission is due to reverse shocks 
crossing the shells. However, this model requires that the ex- 
ternal shock emission is basically turned off. This may need 
conditions difficult to meet. Other authors argue that the jet 
breaks are actually hidden in the optical lightcurves (Curran 
et al. 2007) and/or less evident than expected (Panaitescu et 
al., 2007a, Liang et al. 2007). In Panaitescu (2007b), the au- 
thor proposes a complex scenario, in which the plateau, the 
flares and the chromatic breaks seen in the X-ray lightcurve 
are caused by scattering of the forward-shock synchrotron 
emission by a relativistic outflow, located behind the lead- 
ing blast- wave. Efforts have also been made to reconcile 
the chromatic breaks with the scenario of an unique out- 
flow (Panaitescu et al. 2006a), hypothesizing an evolution 
of the microphysical parameters, including the fractions of 
blast wave energy given to electrons and to the magnetic 
fleld. However, as the authors themselves pointed out, the 
required evolution is assumed "ad hoc", and still lacks a 
self-consistent physical explanation. 

Recently, Oates et al. (2007) have investigated the case 
of Swift GRB050802, one of the bursts in the dataset of 
Panaitescu (2006a), which shows a very evident chromatic 
break. They found that the observed late SED cannot be 
reproduced by models based on single component outflow, 
and proposed a model based on two outflows: a narrow one 
responsible for the X-ray emission, and a wider one that 
powers the optical emission. Both outflows receive continu- 
ous energy by means of shells emitted at late times or in the 
form of Poynting flux. The break in the X-ray lightcurve, 
in this scenario, is interpreted as a jet break, and there is 
no discontinuation of energy injection. The "normal" de- 



cay phase is then a post jet break phase with a slope less 
steep than usual because of the energy injection. The fact 
that the optical lightcurve does not show a break within 
the time of the follow-up observations is naturally explained 
by the lower degree of coUimation of the outflow responsi- 
ble for it. In this paper, we conduct a detailed analysis of a 
sample of other GRBs that are reported to have chromatic 
breaks, showing that this model can potentially interpret 
the observed behaviour. We also discuss how this scenario 
may change our interpretation of the canonical lightcurve 
of GRBs and the deep implications that this change of per- 
spective may have on our understanding of GRB physics. 
This paper is organized as follows. In §[2] we introduce the 
dataset and the data analysis, while in § [3] we present the 
application of the model to the GRB sample. Discussion and 
conclusions follow in § |4] and § (5] respectively. 



2 DATA REDUCTION AND ANALYSIS 

In this work, we reexamine all Swift GRBs with chro- 
matic breaks contained in the sample of Panaitescu et 
al. (2006a), namely GRB050319, GRB050401, GRB050607, 
GRB050713, GRB050802, GRB050922c, in the light of the 
results found by Oates et al. (2007) on GRB050802. We also 
include in our analysis Swift GRB060605, which is another 
example of a burst with a chromatic break and good quality 
data. 

As we will discuss later on, while the X-ray analy- 
sis alone can indicate that our model is compatible with 
the observations, the presence of a second outflow can be 
robustly confirmed only by a joint analysis of the X-ray 
and optical data. In this respect, we note that two bursts 
in the Panaitescu's dataset, GRB050607 and GRB050713, 
have poorly sampled optical data, while for a further one, 
GRB050401, no UVOT data are available because of the 
presence of a bright star in the field of view. For these events, 
we will only consider the X-ray emission, to show that our 
scenario is fully consistent with the observations. 

Once a GRB has been detected by the BAT, Swift im- 
mediately slews, allowing the XRT and UVOT to provide 
prompt simultaneous multi-band data. In the following, we 
describe how XRT and UVOT data are reduced and anal- 
ysed. 

2.1 XRT data reduction 

To determine the X-ray properties of the GRBs, we first 
re-ran the processing pipeline version 2.72 of the Swift 
software. We generated light curves using the software of 
Evans et al. (2007) which supplies the Swift XRT Ught 
curve repositorj{3, and modelled them with a sequence 
of connected powerlaw decays, using minimization. 
In this way we identified the segments of the lightcurves 
corresponding to the lightcurve segments of the canonical 
X-ray lightcurve. We then extracted spectra and effective 
area files (ARFs) for the plateau and post-plateau phases. 
Where the source was piled up, we fitted the source PSF 
profile with SwtfVs known PSF (Moretti et al. 2006) to 

^ http:/ /www. swift. ac.uk/xrt_curves 



Jet breaks at the end of ... 3 



determine the radius within which pile-up is important, 
and used an annular extraction region so that data from 
the piled-up part of the PSF was excluded. If the source 
was not piled up, we used a circular extraction region of 
20 pixel radius (or smaller for faint sources, to maximise 
the signal-to-noise). In some cases, a single light-curve 
segment could cover several decades of count-rate, with pile 
up only being a problem at the start of the segment. In 
these cases we extracted two event lists, using an annular 
source region when pile-up occurred and a circular one 
at all other times, and created separate ARF files for the 
two extraction regions. The event lists were then combined 
using XSELECT and a single spectrum was generated from 
the extracted events; the ARFs were merged using the 
ADDARF tool, and weighting the component ARFs by 
source count-rate. Background spectra were always ex- 
tracted from an annulus centred on the source; these annuli 
were searched for sources, and any found were excluded 
from the extraction region. Where a light curve segment 
spanned multiple Swift observations, separate event lists 
and ARFs had to be generated for each observation; these 
were also combined as just described. Where a spectrum 
corresponding to a specific time was required to produce 
a combined UVOT-I-XRT spectral energy distribution, we 
first determined the count rate C at the epoch of interest 
from the best fit parameters of the light curve, then we 
modified the exposure time in the spectral file so that the 
resulting count rate was equal to C. 



2.2 UVOT data reduction 

UVOT observes the GRB field through a number of pre- 
planned exposures. The automatic target(AT) sequence be- 
gins with a short settling exposure followed by either one 
or two finding charts. UVOT performs observations either 
in event mode, where the position and arrival time of each 
photon is recorded; or, in image mode, where an image is 
accumulated over a fixed period of time. The GRB is ex- 
pected to vary over the shortest timescales during the first 
few hundred seconds after the trigger; therefore, the settling 
exposure and finding charts are observed in event mode. 
The rest of the AT sequence contains a series of exposures, 
in the 7 filters, lasting from as little as 10s through to a few 
thousand seconds. These are observed through a combina- 
tion of event (until ~850s after the trigger) and image mode 
observations. 

The aspect and astrometry for each photon, in the 
case of the event data, was refined following the method 
of Gates et al. (in prep.). The images were processed by 
the pipeline at the Swift Science Data Center (SDG). Any 
images not aspect corrected during the pipeline processing 
were corrected using bespoke aspect correction software. To 
produce lightcurves, the source counts were extracted in an 
aperture which was sized according to the count rate. For 
count rates higher than 0.5 counts per second, a 5" radius 
circle was used, and for count rates lower than 0.5 counts 
per second the source count rates were obtained using a 
3" radius circle, and were then corrected to 5" using the 
PSFs recorded in the calibration files (Poole et al. 2007). 
The background count rates were determined using a circle 
of radius 20", positioned over a blank area of sky near the 



source position. The software used to extract the count 
rates can be found in the software release, Headas 6.3.2 
and version 20071106(UVOT) of the calibration files. In 
order to produce a single optical light curve for each GRB 
in the sample, the lightcurves in each UVOT filter were 
renormalised to that in the V filter. The normalisations 
were determined by performing a simultaneous power law 
fit, in which the lightcurves in the different filters have 
the same slope but were allowed different normalisations, 
in periods in which the lightcurve can be described as a 
powerlaw decay. The count rates from each filter were then 
binned by taking the weighted average in time bins of 
(5T/T = 0.2. 

In order to understand the properties of GRBs of our 
sample, we built the Spectral Energy Distributions (SEDs) 
at two epochs, before and after the end of the plateau in 
the X-ray lightcurve. As for the optical, we used the best fit 
normalisation for each filter to get the corresponding count 
rate at the epoch of interest, by using the best fit decay 
index. The uvottools "uvot2pha" and "ftedit" were used to 
create the spectral files and convert the count rate to the 
value determined in the lightcurves fitting described above. 

2.3 Spectral modelling 

All spectra were fitted in XSPEC 12.3. The X-ray spectra 
were binned to contain a minimum of 15 counts per bin (20 
counts for the brightest spectra), and we used the version 
10 response files (Godet et al. 2008). Some of the plateau- 
phase data comprised both Windowed Timing (WT) and 
Photon Counting (PC) data, in which case the two modes' 
spectra were fitted together with the same model, but a 
(free) constant factor applied to the normalisation. 

Theoretical predictions and observational findings indi- 
cate that the spectral shape of a GRB afterglow is typically 
either an unbroken or a broken powerlaw throughout the 
X-ray and optical bands. The break frequency is the syn- 
chrotron cooling break, uc, in which case the difference in 
the spectral slopes of the broken powerlaw is 0.5. Therefore, 
we jointly fitted the optical and X-ray SED with two models. 
One model consisted of unbroken powerlaw, two absorbers 
and two dust models (zdust in Xspec). The column density 
of one of the two absorbers was fixed to the Galactic value at 
the coordinates of the GRB, given by Kalberla et al. (2005), 
while the value of reddening in one of the zdust model was 
frozen to the value derived from the absorber value, accord- 
ing to relation between Eb-v and the hydrogen column den- 
sity (Bohlin et al. 1978) . The redshift of the other absorber 
and zdust component was fixed to the corresponding burst 
redshift 0. The second model was different only in substitut- 
ing the the powerlaw with with a broken powerlaw, with the 
second spectral slope bound to be higher than the first by 
0.5. In the process of spectral analysis, we tried the Galac- 
tic, Large Magellanic Cloud and Small Magellanic cloud ex- 
tinction laws (Gal, LMC and SMC henceforth). However, 
since in all cases (apart from GRB050802, see below) it has 
been impossible to disentangle among these three extinction 

^ In this regard, all the bursts for which we built the optical and 
X-ray SEDs have their redshift known by spectroscopy 
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laws, in the following we report results obtained adopting 
the SMC extinction law, which provides acceptable results 
in the fits of the extinction laws of the GRB host galaxies 
(Stratta et al. 2004, Schady et al. 2007). For spectral mod- 
elhng of those bursts which only have X-ray data, the model 
was reduced to a single powerlaw and the two photoelectric 
absorbers. 

In the following sections of this paper, we use the con- 
vention F ~ t~°'u~'^ and errors are indicated at 68% confi- 
dence level (c.L). The subscripts "O" and "X" refer to opti- 
cal and X-ray respectively. We will add the labels "1", "2", 
etc to attribute the decay and spectral slope to the relative 
portion of the canonical X-ray lightcurve. The segments of 
the X-ray canonical lightcurve which will thus be ax,i, ax,2, 
ax, 3- The time when the breaks in the X-ray lightcurve oc- 
cur will be indicated as tx,i and tx,2- ax,2 will always be 
the decay slope of the slow decaying segment, ao will be 
the slope of the optical lightcurve. If any break is detected 
in this band, we will define ao,i and ao,2 as the pre and post 
break slope, respectively. The labels Px.i, Px,2 and (ix,-i will 
indicate the spectral energy slopes of the X-ray data only. 
As for the analysis of the SED, in the case of a fit with single 
powerlaw, Pox is the energy index of the spectrum. In the 
case of a fit with broken powerlaw, we shall use two energy 
indeces, which will be referred to as Pox and Pox +0.5 (we 
remind that the difference between the two indeces is fixed 
to be 0.5). Additional "E" and "L" labels indicate if the fit 
was performed before or after the break in the X-ray. 
The results of the temporal analysis of the GRBs are given 
in Table [1] and those of the spectral analysis are reported 
in Tables [31 S] and [S] The formulae we shall be using are 
recollected in Tab. |2l 



3 RESULTS OF GRB DATA ANALYSIS. 
3.1 GRB050319 

The X-ray lightcurve of GRB050319 (Fig [T] top panel) 
shows the typical canonical behaviour, and can be ade- 
quately fitted by a double broken powerlaw model, which 
yields x ~ 113.3 for 113 d.o.f.. The best fit parameters 
are decay indices of ax,i ~ lAS'to'li, ctx,2 = 0.48 ± 0.03, 
ax,s = 1.4lj;g Q^ and break times among these segments of 
tx,i = 300t.ll s and tx,2 = '29.9tli ks. The early, rela- 
tively steep decay is likely the tail of the prompt emission, 
a mechanism that does not involve the forward shock; we 
will therefore ignore this part of the emission hereafter. The 
initial flat decay phase, between tx,i and tx,2, has a spec- 
trum with a powerlaw index Px,2 ~ 1.00 ± 0.03. After tx,2, 
the X-ray spectrum shows marginal indication of softening, 
since the best flt index is Px,3 ~ 1.12 ± 0.07. 

The fit of the optical lightcurve (Fig[T] top panel) with a 
single powerlaw provides a marginally acceptable fit, yield- 
ing = 54.6 with 25 d.o.f. The best fit decay index is 
ao = 0.62 ± 0.02. A fit with a broken powerlaw is slightly 
better, yielding = 43.2 for 23 d.o.f. The F-test indicates 
that the probability of chance improvement is very marginal, 
6.5%. As for the broken powerlaw model, the values of the 
best fit parameters, other than the first slope, are not well 
constrained, if we leave all of them free to vary. We then 
fixed the value of the second slope, forcing it to differ from 



the first decay slope as much ax,2 differs from ax, 3. We 
thus obtained ao.i = 0.58 ± 0.03, ao,2 = 1.52 for the two 
decay indices, and a break time to = 164.4t^g''g^ ks. To find 
a strong upper limit on to, we varied its value while fitting 
the other parameters, until we obtained Ax^ of 9. We found 
that we have to > 51.5 ks at 3 sigma confidence level. There- 
fore, we note that a break in the optical band, if any, takes 
place much later than the break in the X-ray. Our result 
are consistent with those of Panaitescu et al. 2006, in which 
the authors do not find any steepening of the optical band 
emission up to ~ 400ks after the trigger. All these findings 
indicate that GRB050319 has got a genuine chromatic break 
in the X-ray band only at about 30ks after the trigger. 

The SEDs of GRB050319 were built at 20ks and 70ks 
after the trigger (Fig. [U top panel); results of the fit are 
shown in Table [3] For both SEDs, the fit with a cooling 
break in the spectrum yields a better x^ than the fit with 
a single powerlaw, which is nevertheless still acceptable. In 
the following we will discuss both the cases of unbroken and 
broken powerlaws. 

Let us first consider a scenario in which the X-ray and 
optical bands lie on the same spectral segment at 20ks, be- 
low the cooling frequency. This corresponds to the spectral 
fit with a single powerlaw of slope Pox,e ~ 0.84 ± 0.05. 
In this scenario, one should expect that the fiuxes of both 
bands decay with the same slope. We find that the X-ray 
slope observed at early times, ax,2 = 0.48 ± 0.03, is con- 
sistent within ~ 2.4(T with ao,i ~ 0.58 ± 0.03. The aver- 
age decay index of X-ray and optical is a = 0.53 ± 0.02. 
Such a shallow optical decay requires that energy injection 
takes place. The value of the energy injection parameter q 
is linked to the values of the spectral and decay indices {P 
and a) through the expression collected in Tab. [5] (Zhang et 
al. 2006, Panaitescu et al. 2006b). In the case at hand, we 
have q — 0.50 ± 0.06 in the standard hypothesis of a con- 
stant density environment (ISM). The break in the X-ray 
lightcurve at 30 ks is generally interpreted as the cessation of 
energy injection. However, if this were the right scenario, the 
optical emission decay slope would simultaneously increase 
up to Q = 3Pox,e/2 — 1.26 ± 0.08, similar to the X-ray de- 
cay slope. This prediction is not consistent with our analysis. 
Alternatively, if the 30 ks break in the X-ray band were due 
to the transit of the cooling frequency below the X-ray band 
and not to the end of energy injection, the expected post- 
break decay index would be (Tab[2} a = 0.95 ±0.08 whereas 
the observed value is ax,3 ~ lAlt^'Q^- Another possibility 
would be that the cooling frequency is already between the 
optical and the X-ray bands at the time of the first SED. 
This corresponds to the broken powerlaw fits, where we find 
a low energy spectral slope Pox,e = 0.49 ±0.05 at 20 ks and 
Pox,L = 0.581° 'ijj at 70ks. The corresponding high energy 
spectral slopes are set to be higher by 0.5. If the cooling 
break is between the two bands, the only scenario that can 
explain why the X-ray fiux decays slower than the optical, 
before the break at 30ks, is one in which the density profile 
of the circumburst medium is typical of a wind ejected by 
a massive star (with density decreasing as where r is 
the distance from the centre of the explosion and 5 ~ 2; see 
Chevalier & Li 2000). However, even this scenario cannot 
explain the decay slopes of X-ray flux and optical emissions 
after the 30ks break. In fact, the conventional interpreta- 
tion of the canonical X-ray lightcurve is that after the 30 ks 
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break the ejecta do not undergo any further increase of their 
kinetic energy. Without energy injection, in a wind environ- 
ment, the decay slope above the coohng frequency would be 
less steep than that of the optical by 0.25, which is obvi- 
ously not in agreement with our observations. For example, 
the optical slope we would expect is a = |/3 -|- 1/2, where /3 
is the spectral slope in this band. Taking /3 as the weighted 
average of the low energy spectral slopes, the optical decay 
should be ao ~ 1.25 ± 0.08, and the X-ray decay should be 
ctx ~ ao —0.25 — 1.00 ±0.08, which is evidently in contrast 
with our findings. 

In summary, this shows that the steep late X-ray decay 
is not explained if we assume that X-ray and optical are 
originated by the same component. 

We can now demonstrate that the late X-ray break can be 
easily explained as a jet break, under the assumptions that 
the outflow responsible for the X-ray is different from that 
producing the optical emission, and the energy injection rate 
does not change till the end of the observations. Here and 
in the following, we will only consider the simple case of 
side-spreading jet and constant density medium with the 
addition of energy injection (Panaitescu et al. 2006b, but 
see § |4] for a discussion). In such a model, the energy is 
assumed to increase as a simple powerlaw, E oc t^^~''\ and 
the energy injection parameter q does not change with time. 
To compute the value of q, we need the decay and spectral 
slopes, as in previous cases. In the X-ray, decay index is 
ctx.2 ~ 0.48 ± 0.03, while for the spectral index we can take 
the weighted average of the energy index found by the X-ray 
data analysis throughout the whole lightcurve, /3_y = 1.02 ± 
0.03. With these values of parameters and in the case of the 
X-ray band above i^c, we derive (Tab. [2]) g = 0.46 ± 0.06. If 
there were not such energy injection, the decay slope after 
the jet break would become a = 2(5 — 2.04 ± 0.06; but 
the addition of energy into the blastwave flattens the slope, 
leaving the flux decaying with a — 1.31±0.10, a value which 
is within la from the observed one, ax,3 = 1.411q q7. 

In order to compute the size of the beaming angle, 9, of 
the narrow outflow, we use the expression (Frail et al. 2001): 

e = 0.093 (i^) rad, (1) 

where tj^d is the jet break time in days, Ek,52 is the isotropic 
kinetic energy of the outflow, and n is the density of the 
environment in protons per cubic centimetre. 

As we will discuss later on (see § 4), in order for our 
model to hold the kinetic energy in the outflow responsi- 
ble for the X-ray emission should be of order ~ 10% of the 
whole energy of the ejecta (Scenario B, see par. 4). Fur- 
thermore, both the density n of the environment and the 
efficiency rj of the conversion of kinetic energy into gamma 
rays should be moderately low. We assume n = 5 x 10~^ 
and ri ~ 0.01. In order to derive an estimate of the energy 
produced by this burst, we look at the prompt emission flu- 
ence and spectrum. GRB050319 prompt emission between 
20 and 150 keV was fitted by a single powerlaw spectrum, 
with photon index F = 2.1 and had a fluence of 1.1 x 10~^ 
ergs cm~^ (Cusumano et al. 2006). If we assume that the 
prompt emission spectrum of this GRB is described by the 
Band function, with spectral break below 20 keV and a typ- 
ical low energy photon index 1, we find that this burst emit- 
ted 7 X 10^^ ergs in the 1-10000 keV band, on the basis of 



isotropic emission at redshift z=3.24. Under the previous as- 
sumption on efficiency, density and fraction of total energy 
which goes into the narrow outflow, a jet break at 30 ks is 
compatible with a beaming angle of 9n = 0.015 rad. 

We note that, strictly speaking, in Eq[T] we should have 
taken into account that the energy of the ejecta is increas- 
ing during the afterglow. Nevertheless, considering the weak 
dependence of 9 on Ek.52, the value of 6 we found can be 
considered correct within a factor 2. 



3.2 GRB050802 

In the case of GRB050802, we only briefly summarize the 
results obtained by Gates et al. (2007); the X-ray and opti- 
cal lightcurves are shown in Fig[T] (middle panel). The X-ray 
lightcurve breaks from a decay slope of ax, 2 ~ 0.63 ± 0.03 
to a slope of ax,3 = 1.59 ± 0.03, 5.0 ± 0.3 ks after the trig- 
ger. The optical lightcurve is well fitted by a single power- 
law decay with slope ao = 0.82 ± 0.03; the 3a lower limit 
on any possible break in the optical is t=19ks. Two SEDs 
were built at 500s and 40ks after the trigger (Fig. [T] middle 
panel). In the case of GRB050802, the best fit was provided 
by adopting the Gal model. Therefore, for this burst, the ex- 
tinction was determined by applying this law. By applying 
the extinction determined in the early SED to the late time 
SED, it was determined that the late UV/optical emission 
lies above the extrapolated X-ray spectrum. This indicates 
that the optical emission is not produced by the same out- 
flow that is responsible for the X-ray emission, regardless of 
where the synchrotron peak frequency and cooling frequency 
lie. Instead, the double component scenario described ear- 
lier was found to be consistent with the data if the X-ray 
band lies below the synchrotron cooling frequency vc- In 
this case, with the values of parameters ax,2 ~ 0.63 ± 0.03 
and f3x = 0.88 ± 0.04 we can derive q = 0.51 ± 0.06 If 
the break at 5 ks is interpreted as a jet break, the expected 
post-break slope would be (see again Panaitescu et al., 2006) 
Q = 2/3 -f 1 = 2.76 ± 0.08 in case the decay proceeds without 
further energy injection, and a = 1.83 ± 0.17 in case there 
is no cessation of energy injection, which is consistent with 
the observed value of ax,3 within 2a. Results of the analysis 
are shown in Tables [T] and 2] 



3.3 GRB050922C 

A first inspection of GRB050922c data clearly shows a break 
in the optical and XRT lightcurves (Fig. [21 . In order to 
quantify its significance, we fit the lightcurves with a single 
and a broken powerlaw. The early optical emission shows 
some features superimposed on the powerlaw decay, such 
as an evident bump at ~ 150s after the trigger. There- 
fore, we excluded from the fit UVOT data taken during 
the first 200s after the trigger and, for consistency, we did 
this with the X-ray data as well. In the case of the X-ray 
lightcurve, we found that the fit with an unbroken power- 
law yields = 219 for 119 d.o.f., while a fit with a broken 
powerlaw provides an = 141 for 117 d.o.f.. The Ftest 
(Bevington et al. 1969) indicates that the probability of im- 
provement by chance is less than 7 x lO"^'^. For a broken 
powerlaw, the best fit parameters are; initial decay slope 
ax,2 = 1.10 ± 0.02, break time tx,b.2 = 6.45+o ks, and 
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GRB 



ao 



ax, 2 



*x,2 (ks) 



/3x,3 



050319 
050802 
060605 



0.62 ± 0.02 
0.82 ± 0.03 
0.83 ± 0.04 



0.48 ± 0.03 
0.63 ± 0.03 
0.41 ± 0.03 



1.00 ± 0.03 
0.89 ± 0.04 
1.04 ± 0.07 



29.93+2.55-2.80 
5.0 ±0.3 
7.73 ± 0.38 



1.4ll°;°f 1.12 ±0.07 
1.59 ± 0.03 0.88 ± 0.04 
1.93lo:lJ 1.20 ±0.09 



050401 
050607 
050713A 



0.56 ±0.02 0.99 ±0.02 
OMt^ H 1.04 ±0.14 
0.58 ±0.03 1.27 ±0.04 



4.27 ± 0.52 
16 2+'^-* 

7 r;4+0-87 
' -O^-O.SO 



1.44 ± 0.07 0.95 ± 0.07 
1.33l°;Jf 1.17 ±0.20 
1.21 ±0.03 1.02 ±0.05 



Table 1. Results of the analysis of the bursts with chromatic breaks considered in this paper. Prom left to right; burst name, decay index 
in the optical, X-ray decay slope of the plateau phase, X-ray spectral slope in the plateau phase, X-ray lightcurve break time, X-ray late 
times decay index. X-ray late times spectral slope 







no injection 






injection 








13 


a 


a(/3) 




a 




a(/3) 


ISM and spehrical expansion 


I'm < 1^ < l^c 


^ (0.7) 


2^2^(1.05) 


3/3 


(2p- 


-&) + (p+3)q 


(0.38) 


a = (9 1)±'%^"' 




1 (1-2) 


2^ (1.30) 


3/3-1 


(2p- 


-4,) + (v+2)q 
4 


(0.75) 




ISM and jet expansion 




E^i (0.7) 


P (2.4) 


a = 2^ ± 1 


(2p- 


-3)-Kp-|-3)9 


(1.5) 


^ _ (4/3-l)-|-2(/3-H2)9 




1 (1.2) 


P (2.4) 


a = 2/9 


2(p- 


-\-) + {v+2)q 
3 


(1.67) 


2(9- 1) 2{l+q)D 
3 ' 3 


Wind and spherical expansion 


I'm < < I'c 


^ (0.7) 


^ (1.55) 


" 2 


(2p- 


-2) + (p-|-l)q 


(1.13) 


a = f ± (^+/"' 


V > Vc 


f (1-2) 


^ (1.3) 


3/3-1 


(2p- 


-4) + (p-H2)q 

4 


(0.75) 





Table 2. Table with the relations between the value of decay index a and the spectral slope ,i3 in various afterglow models with the 
inclusion of the cases of energy injection. The case of p < 2 is not included, and the self-absorption effect is not discussed. We do not 
consider the case of observing frequencies below Vm- The convention Fu oc t^°'u~^ is adopted here. The temporal indices with energy 
injection are valid only for q < 1, and they reduce to the standard case (without energy injection: Sari et al. 1998, Chevalier & Li 2000; 
with Energy injection: Zhang et al. 2006, Panaitescu et al 2006b) when q = 1. For q > 1 the expressions with energy injection are no 
longer valid, and the standard model applies. The numerical values quoted in parentheses are for p = 2.4 and q = 0.5. 





Fit at 20 ks 




Fit at 70 ks 




Parameters 


Single powerlaw 


Broken powerlaw 


Single powerlaw 


Broken powerlaw 


i3i 
Eb 

132 

Eb-v 
Nh 

Xv 


0.88l»:'^l 

i3.iil:§ X 10-2 

< 0.56 
129.5/108 


0.49 ± 0.05 

2Q+°'^* 

^■^^-0.14 

0.99 ± 0.05 
4.lt| ° X 10-2 
0.62 ± 0.20 
110/107 


0.84 ±0.05 

8.80l° *? X 10-2 
< 0.23 
34.6/24 


0.581^1^ 
0.28«:«« 

l.08l°;l^ 

4.012-9 X 10-2 

< 0.87 
22.2/23 



Table 3. Best fit values of the GRB050319 SED at 20 ks and 70 ks. Nh is expressed in units of 1022 cm-2, the break energy Eb is 
given in keV, and the local extinction Eb_y is in magnitudes. All upper limits are at 90% c.l.. 
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Fit at 20 ks 




Fit at 70 ks 


Parameters 


Single powerlaw 


Broken powerlaw 


Single powerlaw Broken powerlaw 


/3i 


0.86 ± 0.02 


0.89 ±0.04 


0.99 ±0.02 


Eb 




4^^ X 10-3 




P2 




1.39 ±0.04 




Eb-v 


18 ± 2 X 10-2 


18 ±0.02 X 10-2 


0.18 


Nh 


0.26 ±0.04 


0.29 ±0.04 


0.26 


Xv 


120/104 


119/103 


27/15 



Table 4. Best fit values of the GRB0500802 SED at 0.4-1 ks and 35-55 ks. Nh is expressed in units of 10^^ cm-^, the break energy 
Eb is given in keV, and the local extinction Eb—v is in magnitudes. In the case, the fit of the two SEDs was performed by assuming a 
Galactic extinction law (all results are taken from Oates et al. 2007). 



late decay slope ax, 3 ~ l-iS'^^'g^. For the optical lightcurve, 
a single powerlaw fit of the renormalised V, B and U band 
lightcurves gives ~ 110.2 for 18 d.o.f., whereas a broken 
powerlaw gives ~ 24 for 16 d.o.f.. In the latter case, the 
best fit parameters are ao,i = 0.77 ± 0.03, to,b = 6.23lQ;gg 
ks, and ao,2 = 1-20 ± 0.05. As we can see, from our re- 
analysis and new reduction of the X-ray and optical data, 
the break times in the two bands turn out to be consistent 
with each other within la, suggesting that the break in the 
X-ray lightcurve should not be considered as achromatic, in 
contrast to what was suggested by Panaitescu et al. (2006a). 



3.4 GRB060605 

The Smft GRB060505 also shows a canonical X-ray 
lightcurve, with an initial steep decay, a shallow plateau 
and finally a steep decay (Fig. [1] bottom panel). The decay 
slopes of the three segments and the two break times are 
ax,i = "i.GStoil, tx,i = 164.5l?^:^s, ax,2 = 0.41 ± 0.03, 
tx,2 = 7.73 ± 0.38 ks, ax.a = 1.93+15 {J. There is no 
evident strong X-ray spectral evolution, since the X-ray 
energy index in the plateau and in the steep decay are 
Px,2 = 1.04 ± 0.07 and /3x,3 = 1-20 ± 0.09, consistent 
within la. In the optical, GRB060605 shows a wide peak 
at few hundreds seconds after the trigger, which is likely 
to be the beginning of the forward shock emission (Oates 
et al. in prep.). In fitting the optical lightcurve (Fig. [T] 
bottom panel), we considered all the datapoints taken after 
500s from the trigger. The single powerlaw model provides 
a marginally acceptable fit, with = 28 for 13 d.o.f.. 
We then tried a broken powerlaw model, which gives a 
much better fit with ~ 10-6 for 10 d.o.f.. The value of 

but it is not well 



the late decay slope is 0:0,2 ~ 3.3_ 
constrained; we can infer that it has a lower limit of 1.4 at 
95% C.L.. The best fit values of the other parameters are 
ao,i = 0.85 ± 0.04 and to = '^S-^ttl ks. The 3a lower limit 
on the break time in the optical, calculated as in the case of 
GRB050319, is to = 12.3 ks. Ferrero et al. (2008) present 
a dataset in which the optical afterglow is well detected till 
~ 1 day after the trigger, and their data show an evident 
break occurring 23.3 ks after the trigger, with a late decay 
slope ao,2 = 2.56 ± 0.16. We note that our best fit values 
are consistent with those of Ferrero et al. (2008). Thus, we 
can conclude that a break is present in the optical, but it 



is inconsistent with tx,2. Ferrero et al. (2008) suggest that 
the different break times might be caused by some fiaring 
activity in the X-ray band that occurred around 6ks after 
the trigger. These flares would have led to the conjecture 
of an X-ray afterglow decaying shortly thereafter (see their 
paper for more details). We will rather investigate the 
scenario in which GRB060605 has a genuine chromatic 
break. For this GRB, we built up the SEDs at 5ks and 
20ks; the values of the best fit parameters are reported 
in Table [5] As we can see, we cannot distinguish between 
the single powerlaw and the broken powerlaw spectral fit 
on statistical basis, since both models provide a similar 
reduced x^- However, an unbroken powerlaw model is ruled 
out by the fact that X-ray and optical decay slopes are 
inconsistent at 7a level at 5 ks and IOct level at 20 ks. We 
are thus left with a scenario in which the spectrum is a 
broken powerlaw at 5 ks and 20 ks. Furthermore, we have 
to assume a wind circumburst environment for the same 
reasons quoted for GRB050319. We reiterate that, in this 
environment, the cooling frequency is supposed to increase. 

If we fit the two SEDs with a broken powerlaw and 
restrict the break energy between 0.005 and 1 keV, the low 
energy spectral indeces are f3ox,E = 0.54 ± 0.07, Pox,l = 
0.71 ± 0.09, at 5 and 20 ks respectively. The break energy 
at 5ks is 0.008 keV, with a 1 a positive error of 0.032. This 
break energy value is near the minimum allowed value of 
0.005 keV; we were not able to find an la negative error. 

We note that the low energy spectral indices are con- 
sistent within 2a. We assume an average low index Pox ~ 
0.60 ±0.06 and a high energy index ^^^^ + 0.5 = 1.10 ±0.06 
respectively. The first index has got to be that of the Op- 
tical band. In the usual interpretation of the canonical X- 
ray lightcurve, the break at 7.3ks corresponds to the end 
of energy injection into the ejecta. If this is the right sce- 
nario, in a wind density profile, the optical emission de- 
cay index should be higher than that of the X-ray emis- 
sion by 0.25. For example, we should observe an optical 
decay slope ao = f^ox + 1/2 = 1.40 ± 0.09 after the 
end of the injection; the X-ray fiux decay index ought to 
be ax ~ Qo — 0.25 — 1.15. These predictions are clearly 
inconsistent with the observed behaviour. The X-ray flux 
would decay faster than 1.15 if the cooling frequency moved 
above the X-ray band, but in such a case the X-ray decay 
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slope would be consistent with that of the optical, which is 
inconsistent with observations, as stated above. 

We try now to apply our model to interpret the be- 
haviour of the X-ray emission for this burst. Again, the idea 
is that the plateau, extending till 7.3 ks after the trigger, 
is due to forward shock emission of ejecta expanding like 
they were spherical, with the contribution of energy injec- 
tion. For this burst, we suppose that the X-ray band remains 
below the cooling frequency. In fact, by using ax.2 and the 
weighted average energy index = 1.10±0.06, from Tab. [2] 
we derive q = 0.20 ± 0.06. Assuming that the end of the 
plateau phase is due to a jet break with side expansion, the 
predicted decay slope post break would be a = 3.20 ± 0.12 
or a = 1.48 ±0.20 in case of cessation or continuation of the 
energy injection, respectively (see Tab[2}. Again, the sec- 
ond value is consistent with the observed result at 2a level 
( Ox, 3 = 1.931q5^J). -fri order to compute the opening angle 
On of the outflow responsible for the X-ray emission, we can 
follow the same procedure as GRB050319 after estimating 
the total emitted energy. According to Sato et al.(2006), the 
fluence in the 15-150 keV band of GRB060605 is 4.6 x 10"^ 
erg cm~^, while the spectrum is best fitted by a simple pow- 
erlaw with photon index Fi ~ 1.34. Since this value suggests 
a high energy spectrum below the break energy (Band et al. 
1993), we can assume that the break energy is occurring just 
above the BAT bandpass. Assuming that the high energy 
photon index is F2 = 2.3 (the average value for this param- 
eter foUowing Band et al. 1993) and redshift z = 3.8, we 
find that the isotropic energy emitted between 1 and 10000 
keV is _E ~ 3.5 x 10^^ ergs. The next step is to estimate 
the kinetic energy of the ejecta and which fraction of it goes 
into the narrow outflow. Now, in the case of GRB060605, 
a possible jet break occurs in the optical not much later 
than the jet break in the X-ray, Eq. [T] indicates that the 
opening angle 6w of the outflow responsible for the opti- 
cal emission and 9n could be close. Now, in our modelling 
(see section 4 for details), it is intrinsically assumed that 
we have emissions from spherical portions of two outflows, 
and the emitting surface of the narrow outflow, responsible 
for the X-ray emission, is much less than the surface of the 
wide outflow, which is producing the optical emission. The 
approximation can hold if the beaming angles of the two 
outflows are different enough. A way we can reconcile our 
interpretation with the features of GRB060605 is by assum- 
ing that the energy in the narrow component En is much 
higher than the energy Ew carried by the wide outflow. In 
our theoretical discussion, we have found that solutions with 
En — SOEw are possible (Scenario A", see section 4). This 
solutions applies in cases of density n ~ 1, and efficiency of 
conversion of kinetic energy of the ejecta into 7-ray emission 
rj — 0.2. We thus derive that the kinetic energy of the ejecta 
is ~ 1.8 X 10^^ ergs. Now, if we apply this ratio of energies 
and this density to GRB060605, then we derive, by using 
Eq. [T] that the narrow outflow should have an opening an- 
gle 9n = 0.019 rad. The outflow responsible for the optical 
emission should have 9w ~ 0.05. 

3.5 GRBs with X-ray data analysis only 

All the bursts for which we built the optical and X-ray SEDs 
have their redshift known by spectroscopy, while the follow- 
ing other objects in our sample do not have known redshifts 



(except 050401). However, since they are studied in the X- 
ray band only, the lack of a redshift basically does not affect 
our results and conclusions. 

GRB050401 - A break is evident in the X-ray lightcurve 
of this GRB: the decay slope changes from ax, 2 = 0.56±0.02 
to ax,3 = 1.44 ± 0.07 at tx,2 = 4.27 ± 0.52 ks. There is 
not strong spectral evolution throughout the whole obser- 
vation, since the spectral index is always consistent with 
fHx ~ 0.99 ±0.02. Again, if the X-ray band is below i/c, then 
the energy injection parameter would be g = 0.39 ± 0.03. If 
the outflow responsible for the X-ray emission underwent a 
jet break without energy injection, the predicted slope of the 
flux decay would he a — 2.98 ± 0.04, which is inconsistent 
with the value we observe. However, in the presence of en- 
ergy injection the predicted value is a = 1.74 ± 0.09, which 
is consistent with the observed X-ray decay slope at ~ 2.3a 
level. In order to compute the beaming angle of the outflow 
responsible for the X-ray emission, we need to make some 
assumptions. We will assume that the Energy of narrow out- 
flow responsible for the X-ray emission is 10% that of the a 
wider outflow that produces the optical emission (Scenario 
B), and an efficiency r] — 0.01 and a density n = 5 x 10""^. 
We have E-y = 3.5 x 10^^ ergs (Golenetskii et al. 2005). 
With these assumptions for density, efficiency and ratios of 
kinetic energies the jet beaming angle of the narrow compo- 
nent turns out to be 9n = 0.006 rad (Eq. [T}. 

GRB050607 - This burst exhibits an evident break in 
the X-ray lightcurve, since its decay slopes change from 
ax,2 = 0-54:to°l to ax,3 = l-SStolU at 16.2^^;^ ks. The 
X-ray spectrum does not show evidence of evolution at the 
break time and has an average energy index of /3x = 1.07 ± 
0.11. Assuming that the X-ray band is above the cooling 
frequency, the values of f3x and ax, 2 imply q — 0.59 ± 0.23 
(Tab. [2}. Without late time energy injection, the subsequent 
jet decay slope would be a = 2.14 ± 0.22, while with energy 
injection the predicted value is a = 1.57 ± 0.41. The latter 
is consistent with the observed value of Qx,3, within ~ la. 
In order to derive the beaming angle of the narrow outflow, 
we need an estimate of the burst energetics. Since the red- 
shift of this burst is presently unknown, we adopted z = 2.5 
(i.e. about the average Swift GRB redshift, Jakobsson et al. 
2006) and a prompt emission spectral index estimated by the 
Band function, with a high energy photon index of ~ 2.1 in 
the energy band from 15 keV to 10000 keV and of ~1 below 
15 keV (Pagani et al. 2006 report 1.83±0.14 in the range 15- 
150 keV). Under this hypothesis, the energy emitted by the 
burst would be ^ 2.4 x 10^^ ergs. We can assume that 90% 
of the kinetic energy of the outflows is carried by the broad 
one, and we can take Vc below the X-ray band (scenario B, 
section 4); other assumptions are n = 5 x 10"'', r\ — 0.01. 
With these hypothesis in place, we obtain a beaming angle 
of 9n = 0.023 rad. 

GRB050713A - the X-ray lightcurve of this burst shows 
a break at tx,2 = 7.54lQ gQ ks, after which the decay slope 
increases from ax, 2 = 0.58 ± 0.03 to ax, 3 = 1.21 ± 0.03. 
The spectral index, throughout the whole observation, is 
/3x = 1.17±0.03. The energy injection parameter, again for 
the case of X-ray band above i/c, is g = 0.38 ± 0.06. The 
expected slope at late times would be a = 2.34 ± 0.06 or 
a = 1.44 ±0.11 in case of cessation or continuation of the 
energy injection process, respectively. The latter is consis- 
tent with the observed decay slope in the X-ray band within 
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Fit at 5 ks 




Fit at 20 ks 




Parameters 


Single Pow. 


Broken Pow. 


Single Pow. 


Broken Pow. 


/9i 

Eb 

132 

Eb-v 
Nh 

Xv 


1 01 +°-''^ 

7.64l2-i« X 10-2 
< 0.96 
42.0/42 


0.54 ±0.07 

0.008+0032 

1.54 ±0.07 
0.008+0-032 
49+0-24 
41.4/41 


1.16 ±0.09 

12.5+0-5 X 10-2 
1 90+0.46 

28.4/27 


0.71 ±0.09 

„ o~r:+0.215 
U.O(0_Q 

1.21 ±0.09 

< 0.122 

1 48+0-'^! 

26.4/26 



Table 5. Best fit values of the GRB060605 BED at 5 ks and 20 ks. Nb is expressed in units of 10 cm , the break energy Eb is given 
in keV, and the local reddening Eb-v is in magnitudes. Upper limits on column density and reddening are at 90% c.l.. 



GRB Observed ax,3 Predicted ax,3 



050319 

050802 
060605 



050401 
050607 
050713A 



1+0.08 
-0.07 

1.59 ±0.03 
,+0.11 
-0.10 



1.41"' 



1.93^ 



1.31 ±0.10 
1.83 ±0.17 
1.48 ±0.20 



1.44 ±0.07 
1 oq+0.16 

^••^"^-0.11 

1.21 ±0.03 



1.74 ±0.09 
1.57 ±0.41 
1.44 ±0.11 



0.015 

0.017 
0.019 



0.006 
0.023 
0.008 



Table 6. GRBs with chromatic breaks considered in this work. The table shows the late decay slope observed in the X-ray, the slope 
predicted by our model, and the inferred values of the beaming angle for the narrow outflow. In the case of GRB050802, values axe taken 
by Gates et al. (2007). 



2(T. To calculate the beaming angle of the narrow outflow, we 
made again an assumption on the (currently unknown) burst 
redshift. By using z = 2.5, and taking the values of fluence 
and spectral parameters published in Morris et al. (2007), 
we infer an isotropic 7-ray energy of ~ 2.5 x 10"'''^ ergs. 
With the same assumptions made for GRB050607, we ob- 
tain On = 0.008 rad. 



ent degrees of collimation. The narrow outflow generates 
the X-ray emission, while the wide one the optical. Both 
emissions are due to the usual forward shock, which has a 
synchrotron spectrum consisting of powerlaws connected at 
particular frequencies (Sari, Piran & Narayau 1998), i.e. the 
synchrotron frequency and the cooling frequency i>c- In 
this paper, we use the expressions of Vm, I'c and of the peak 
flux Fmax as determined in Zhang et al. (2007) for a constant 
density medium: 



4 DISCUSSION 

Results reported in the previous section show that a single 
outflow model cannot explain the behaviour of the GRBs 
with chromatic breaks we have considered. Instead, we found 
that if the X-ray flux is attributed to ejecta which are de- 
coupled from those responsible for the optical, the observed 
behaviours of these GRBs can be explained. In the theo- 
retical modelling of GRBs, a double component outflow has 
already been put forward, even before the launch of Swift 
(e.g. Bcrgcr et al. 2003, Peng et al. 2005). It has been invoked 
to explain the complex temporal behaviour of X-ray and op- 
tical emissions of the exceptional GRB080319B (Racusin et 
al. 2008). In this section, we would like to explore the via- 
bility of the two-component jet model with the important 
addition of a continuous energy injection, from a theoretical 
point of view. 

The basic picture is based on ejecta with two differ- 



1600(1 -I- z)D28e-J^_^EK,52n^'' /uJy 

,12 ' 



3.3 X 10' 



p-1 



il + z) 



1/2 1/2 2 

^B,-2^e,-l 



-3/2 - 



6.3 X 10^''(1 + zy^^'^ts-^^lE^^'n'^t, 



1+-1/2 

d 



Hz, 



(2) 



where z is the redshift, 7328 is the luminosity distance in 
units of 102* cm, eB,-2 and Se.-i are the ratios between 
the magnetic/electron and kinetic energy of the ejecta (in 
units of 10~^ and 10~^ respectively), i?52 is the isotropic 
kinetic energy as measured in the observer rest frame and 
normalized to 10^2 ergs, n is the particle density in cm~3^ 
p is the index of the the powerlaw energy distribution of 
radiating electrons, and ta is the observer time in days. 
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By taking z — 2.5 (as for an average Swift GRB, see 
previous sections) and a cosmology with Ho = 71, — 
0.3, A = 0.7, gives D28 ~ 6.2. We adopt a typical value 
of J3 = 2.4, which gives an energy index between Vm and Vc 
of l3 = {p~ l)/2 — 0.7. Below Vrn we assume a standard syn- 
chrotron spectrum rising with /3 = —1/3. In order to take 
into account the energy injection, we assume that the lumi- 
nosity of the GRB central engine scales as L cx: with 
a typical value of g = 0.5 (Zhang et al. 2006). This corre- 
sponds to an increase of kinetic energy of the ejecta of the 
kind E cx t^^~'''' = All these assumptions allow us to 
recalculate the coefficient in the formulae of [5] and change 
the time dependencies, taking into account the increase in 
energy. We obtain 



Fma:, = 2.55 X 10^ e]^^_^E52,on ^h]/\jy 

12 1/2 2 ^p,l/2,-5/4. 



^™ = 2.1 X We;/;_2<_i£;5'4t)*d Hz 

A A inl4 -3/2 jp-1/2 -1 -3/4-u- 
Uc = 4.4 X 10 £5,^2-2^52,0 *d Hz 



(3) 



where £52,0 is the isotropic kinetic energy at 300s after the 
trigger. We chose this time because it is typically from 300s 
that the slow decline phase is observed in both the X-ray 
and optical afterglows. Furthermore, we require our scenario 
to work up to 0.1 days after trigger, since it is typically 
around ~ 0.1 days that the plateau phase ends. To distin- 
guish between the narrow and wide component, we use the 
pedices "n" and "w" respectively, while "O" and "X" indi- 
cate the optical and X-ray band. For the optical and X-ray 
frequencies, we used the values uo ~ 5.5 x 10^** Hz and 
r^x = 10^* Hz, respectively. Therefore, for instance, fo,w is 
the optical flux due to the wide component. 
In the following treatment, we shall be discussing six pos- 
sible scenarios. In order for our model to work, the nar- 
row/wide component should not contribute significantly to 
the optical/X-ray flux. We translate this "condition of non- 
interference" by requiring that the optical flux of the narrow 
component is at maximum one half of that of the wide one, 
and a similar condition for the X-ray band. The six difi'erent 
scenarios we are considering refiect six different possible hi- 
erarchies between the various frequencies. Scenarios A and 
B deal with the case in which both and Ucm lie above or 
below the X-ray band, respectively. The next two cases. A' 
and B', are a variant of the previous ones, in which Uc,m and 
Ucn do not lie on the same side with respect to the X-ray 
frequency. Cases A" and B" show the same arrangements 
of frequencies as A' and B', but the synchrotron peak fre- 
quency of the narrow component is below the optical since 
the beginning of observations. Our data do not allow to dis- 
tinguish between the cases A, A', and A" (or B, B' and B"). 
We require i'm,w < vo and Um,n < ^x, consistently with 
the absence of an increase in the optical and X-ray flux at 
early times in the datasets we have analyzed. All scenarios 
are summarized in Fig. |3] In § 14.71 we discuss the extension 
of validity of the conditions we pose after 0.1 d. 



4.1 Scenario A 

The conditions to apply in scenario A are: 
fo,w > 2f„,o at 0.1 d after the trigger. 



fx,w ^ 2 f^^^ 

i^m,w < 1^0 at 300 s after the trigger. 
i^m.n > 1^0 at 0.1 d after the trigger. 
i^x < Vc.w at 0.1 d after the trigger. 
i^x < Vc.n at 0.1 d after the trigger. 
Vm,n < i^x at 300 s after the trigger. 



(5) 

(6) 
(7) 
(8) 
(9) 
(10) 



It is easy to verify that, if the above conditions are satisfied 
at the time indicated, they are also valid for the whole in- 
terval in which we are interested, i.e. between 300s and 0.1 
days after the trigger. Since the second condition describes 
the evolution of the flux below the cooling frequencies for 
both components, time dependencies cancel out. The first 
condition (Eq.|3)) can be written as 



^1/2 p 



5.5 X 10 



14 



9 1 V inl2,i/^ ,2 pl/2 |-|i_5/4 



> 2<;s,l2,n'^52,0,n 



.1/2 .2 

5.5 X 10" 



9 1 V 1012.1/2 2 pl/2 n-,-5/4 I 

which, after some iterations, can be rearranged into 

0.85 Tpl.35 1.4 ^ on c; 1/3 Ei5/6 -2/3 /^nx 
eS,-2,m-C'52,0,ii.<;e,-l,tu > ^a-O'^B , _2,n -'^52,0,n ^^e, - 1 ,n • U^i 

Similarly, the second condition (Eq.[5l) can be expressed 



1/2 



,,-^52,0, 



1 



< 2^S.-2.n£^52,0,n , ^ ^ 

which simplifies to 

0.85 r,1.35 1.4 , 1 0.85 T-,\.-ih 1.4 

£B,-2,m-C'52,0,m£e,-l,tu < g ^-B ■ -2,n-C'0,52 

From these two inequalities we have 




(1 



(14) 



(15) 



(4) 



E52,0,n > 2.73 X 10 ^B,-2.n^e,-l,n ■ 

We can now obtain a constraint on E52,o,n from Eq. [9] 
E52,o.n < 6.1 X 10-'^esf_2,„n"^ , (16) 
which, substituted in Eq. 1151 gives 

€B,-2,n < 4.8 X lO'^^Ee.-l.nn"' . (17) 

By substituting Eq. [17] into Eq. [15] and after some ma- 
nipulating, we have 

-B52,o,n > 5.4 X 10^e-?.i_„n. (18) 

From Ea. ll8l we can infer that the value of ee,-i,n must 
be very high, in order to avoid an unreasonable value for the 
energy of the narrow outflow. By assuming ee,-i,n = 3.3, 
i.e. the maximum value (which is given at equipartition), 
we obtain i552,o,n > 4 x lO^n. With this value of ee,-i,n, 
a constraint on es can now be obtained from Eq 1171 by 
assuming n = 0.01 it gives es,-2,n < 5 x 10~^, which is a 
very low value. 
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We try now to derive some constraints on the physical 
parameters of the wide component. By solving Eq.[T2]for the 
parameter Es2,o,n and substituting it into Eg. 1141 we derive 



< 2.3 X 10 ^6B,_2,ne^;!.\,„ . 



(19) 
(20) 
(21) 



52,0,m'=S,-2,tii'^e,-l,tu 

which can be combined with Eq. |S] 

E;2i'',^tt/X^n-^ > 404.4 
to obtain 

IT ^ /I V, in6 0-54 -1.29 -0.47 -3.8 

J^52,0,w > 4 X iO n ^e.-l,w^B,-2,n^e,-l,n ■ 

Under the previous assumption of ee,-i,n = 3.3 and 
taking eB,~2,n = 10"^, we derive £^52,0,™ > 3.5 x W'<^~]ll^^- 
Again, the fraction of the energy given to the electrons must 
be close to equipartition, in order to avoid very high values 
of the energy of the wide component. If ee,-2,in ~ 3.3, we 
obtain £^52,0,™ > 8 x 10^. It is physically implausible to have 
a value of £52,0,™ much higher than this lower limit. Apart 
from these caveats, we can now show that the set of inequal- 
ities assumed within scenario A cannot be simultaneously 
verified. In fact, Eq.|6] reads 



1/2 



-2,m'^<i,-l,wE^2,0.w 



1/2 



< 0.221 . 



(22) 



With the values we just obtained for €e,-i,w and £52,0,™, 
Eq. [22] requires an extremely small value of the magnetic en- 
ergy, tB,-2,w < 10~^. On the other hand, by substituting the 
lower limits on £52,0,^, £52,0,™ and ee,-i,n = ee,-i,™ = 3.3 
into Eq.[T2]and 1141 we can derive the following inequalities 



8.8 X 10^es|*!2,n > 10 X lQ'^e'B%.w > 2 x 10^ e^''!. 



2,n ' 



(23) 



and, in turn, a lower limit on 6b, — 2,71 > 7.7 X 10"^ By 
substituting this value in the right member of Eq. [23] gives 
es,-2,m > 3.3 X 10~^, which is in contradiction with what 
was derived from Eq. 1221 

4.2 Scenario A' 

We now consider a variant of the previous scenario, in which 
the cooling frequency of the wide component lays below 
the X-ray frequency but above the optical band. As already 
mentioned, our model cannot distinguish between this sce- 
nario and that described in the previous subsection. The set 
of conditions expressed in Eqs.|3]-[TD]are modified as follows: 



fo,w > 2/0, n 


at 0.1 d after the trigger. 


(24) 


fx,w < ^fx,n at 300 s after the trigger. 


(25) 


l^m,w < vo 


at 300 s after the trigger. 


(26) 


VO < Vc,n} 


at 0. 1 d after the trigger. 


(27) 




at 300 s after the trigger 


(28) 


^ ^m,n 


at 300 s after the trigger 


(29) 


Vm,n > VO 


at 0. 1 d after the trigger 


(30) 


^x *^ ^c^n 


at 0. 1 d after the trigger . 


(31) 



Notice that Eq. [301 must now hold at 0.1 days after 
the trigger, while Eq. [25] must be satisfied at 300s. These 
conditions are translated into the following inequalities: 

0.85 17,1.35 ,1.4 ^ on c;.l/3 rT.5/6 -2/3 



<^B,-2,wJ^52,0,w<^e,-l.w > '^^■O<^B,~2.n^52.0,n''e.-l,n > 



(32) 



0.1 1.4 j-,1.1 ^ n 0.85 1.4 pl-3^ 

^B,~2,w^e,-l,-w-'^52.0,w < ^•O0£s,-2,n£e,-l,ri-C'52,0,n"' 



1/2 



1/2 



-2,iD^e,-l,™-'^52,0,tu 



E. 



1/2 



< 0.22 



3/2 p-1/2 ,_-l 



■=3,-2,111-^52,0,11) 
-3/2 



> 0.22 



■=3,-2,111-^52,0,11) 



n-' < 32.5 



1/2 



2,n'=e,-l,n-C'52,0,i 



1/2 



-2,n'=e,-l,n-^52,0,n 



E. 



1/2 
52,0, 

1/2 



< 4 X lO'' 
> 14.7 



-3/2 p-1/2 ^-1 



> 4 X 10^ 



"=S,-2,n-^52,0,n' 

By combining Eq. [34] with Eq. [36]we obtain 
£52,0,11) < 0.35ejli ^n. 



(33) 
(34) 

(35) 
(36) 
(37) 
(38) 
(39) 

(40) 



In order not to restrict ourselves to solutions with very low 
£52,0,11), we will assume a quite small value of ee,-i,iD 

Rearranging Eq. 1361 we derive a constraint on es,_2,tD: 



£3,-2,1. > 0.10£ '/^n-^/^ , 



while solving Eq. [38] for es. 
Eq. 1391 we obtain 



£5 



> 1.28 X 10 e: 



(41) 

-2,n and substituting it into 
(42) 



From this last equation we infer that £e,-i,n cannot be 
very far from the maximum value of 3.3, achieved at equipar- 
tition, to avoid extremely high values of energy of the narrow 
outflow. If ee,-i,n = 3.3, then £52,0,11 > lO^n. 

Finally, from Eq. [35]we obtain a constraint on £3,-2,11 



es,-2,n < 0.019£, 



-1/3 -2/3 

o n ^ ^ 



(43) 



Now we will show that the previous set of inequalities 
can be solved simultaneously by assuming not unreasonable 
values of the physical parameters, provided that we limit 
ourselves to a scenario in which the circumburst density is 
relatively small, n ~ 10"'^, which makes the set of conditions 
easier to meet. Let us first assume that ee,-i,n = 3.3 and 
that the energy of the narrow component is, £52,0,11 ~ 5. For 
the chosen values of these two parameters, Eq. 1381 requires 
£3,-2,11 > 0.6. Taking £3,-2,^ = 1, from Eq.[3i]and Eq. [32] 
we can derive a lower limit on the values of £52,0,™- In fact, 
Eq. [35] becomes 



0.85 1.4 17)1.35 , rn 

£3,-2,iD£e,-l,iu-C'52,0,ii> ^ OU , 



which, combined with Eq. 1341 gives 

£3^-2, id£52, 0,11) > 150 . 



(44) 



(45) 



Now, the highest value possible of £3,-2,11; is reached at the 
equipartition, £3,-2,™ = 33; in such a case, £52,0,111 > 25. 
By assuming a more reasonable value of 10 
(Panaitescu & Kumar 2001b, Yost et al. 2003) give instead 
£52,0,11) > 50. 

By rearranging Eq. 1321 and 1341 we derive 



-4 -0.5 ^ Q V 

£e,-l,iu£3,-2,iD > o X iU 



(46) 

-2,11) chosen above, implies 



which, with the value of £3,, 

£e,-l,iD < 0.1. 

Let us assume the following series of parameters: n — 
3 X 10"^ £e,-i,n = 3.3, £3,-2,™ = 10, £<7,-i,id = 0.04 and 
£52,0,111 ~ 600. As we will explain in the following, larger 
values of £52,0,11) are implausible, since they translate into 
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unphysically high values of kinetic energy of the wide com- 
ponent at late times. Moreover, since our model requires 
a substantial difference in the beaming angles of the wide 
and narrow components, then the difference in the respec- 
tive kinetic energies must not be too large. We then assume 
Ez2,o,n ~ 0.20i?52,o,iu = 120, and tB,-2,n = 0.1 (the latter to 
satisfy Eg 143 1). This set of parameter values satisfies all the 
required inequalities. We note that the narrow component 
has "standard" values of the two e's (Panaitescu & Kumar 
2001a, 2001b). The wide component, instead, should have an 
inefficient conversion of shock energy into electron energy 
and a very efficient conversion of shock energy into mag- 
netic field. Furthermore, the wide component should carry 
a high amount of energy, since Ek.w is already as high as 
6x 10^** ergs 300 seconds after the trigger, and it increases, in 
our model, as f^^'^ . As mentioned above, EK,m should not be 
much higher than this value. For example, if the initial value 
of the wide component kinetic energy is _Eo,™ = 10^^ ergs, 
this quantity would become as large as -B ~ 3 x 10^^ ergs 4 
days after the trigger. This very high value would likely pose 
a energy budget problem for the central engine of the GRB. 
If GRB optical lightcurve undergoes a jet break several days 
after the trigger (Frail et al. 2001), for these very high values 
of kinetic energy the beaming correction would be ~ 10"** 
(see Eq. [1]). If, instead, Ea^^ = 6 x 10^'*, the kinetic energy 
of the wide outflow would approach lO^'^ ergs 1 day after the 
trigger, and 2 x 10^® ergs 4 days after the trigger. If corrected 
for the beaming factors seen above, the energy of the wide 
component would be of order of lO^'^ ergs which, although 
high, is still acceptable according to GRB theoretical mod- 
els. The large majority of ejecta kinetic energy is carried 
by the wide outflow. Since, in the prompt emission phase, 
the GRBs emit isotropically around 10^^ ergs in gamma ray, 
values of efficiency rj as low as a fraction of percent should 
be assumed (Zhang & Meszaros 2004), at least for the wide 
outflow. 

A possible limit of scenario A' is that, even by assum- 
ing a value of the circumburst medium density as low as 
~ 3x 10~^, it requires a certain degree of fine tuning between 
the parameters. The inequalities required by this scenario 
can be solved also for slightly larger values, i.e. n ~ 10"'^, 
but the allowed region in the parameters space becomes 
smaller and even finer tuning is needed. 

4.3 Scenario A" 

We will now explore a variant of scenario A', which is ob- 
tained by placing the synchrotron peak frequency of the nar- 
row component below the optical band. This condition must 
now hold since 300 s. The new dataset of inequalities reads 



fo,n, > 2/o,„ (47) 

fx,w<^fx,n at 300 s after the trigger. (48) 

Vm.w < (49) 

I'D < i^c,w at 0.1 d after the trigger. (50) 

> Vc,w at 300 s after the trigger (51) 

i'm,n < i^o at 300 s after the trigger (52) 

i^x < i^c,n at 0.1 d after the trigger . (53) 



We note that the inequality fo,w > '2fo,n now has no 
requirement on time, since it deals with fluxes in the same 
spectral regime. However, its expression will have to change 
from the previous scenario. Eq. [52] also changes. We have 

0.85 i-,1.35 1.4 r, 0.85 rril.SS 1.4 /rA\ 

fS,~2,iD-C'52,0,ii>^e,-l,iu > ■^^B,~2,n-'^52,0,n^fE,-l,n ) [O'i) 

0.1 1.4 ipl.l ^ o 0-85 1.4 1^,1.35 1/2 /rr\ 
^B,-2,w^e,-l,w-l^52,0,w < ^.O0es__2 

4'',-2,™^e.-i.™S52,o,™ < 0.22 (56) 

^~B%l,^E;,]i%n-' > 0.22 (57) 

^s'-iAT^^^' < 32.5 (58) 

4''.-2.n^e,-l.„i?52,0,„ < 0.22 (59) 

^~B%lnE;2Tnn'' > 4 X w\ (60) 

In this scenario, ec._i,n should not be so high as 
in other cases. Equations [JD] and 31] still apply. In this 
scenario, it is possible to have values of £^52,0,™ much 
lower than in the previous scenarios and well below 
E52,o,n- In fact, these values meet all the posed condi- 
tions: E52,0,n = 100, ee,-2,n = 0.25, es,-2,n = 2 X 10"^, 

E52,o,m = 3, ee,-i,w — 0.25, eB,-2,w — 2, n = 0.5. This 
fact has important consequences. In our modelling, it is 
intrinsically assumed that we have emissions from spherical 
portions of two outflows, and the emitting surface of the 
narrow outflow is much less than the surface of the wide 
outflow. This approximation can hold if the beaming angles 
of the two outflows are different enough. If 9w — On 
the emitting surface of the wide outflow would be better 
approximated by a ring rather than a portion of spherical 
surface. This conflguration would lead to a behaviour of 
the optical emission which is different from that described 
in our scenario. Now, in the previous scenarios, any break 
in the optical should be much later than the chromatic 
break in the X-ray, otherwise, from Eq. [T] we would have 
indeed drawn that 9w ~ 6m- This stems from the fact 
that in all previous scenarios £52,0,™ is much higher than 
E52,o,n- However, in Scenario A", it is £52,0,™ — 0.03i552,o,n. 
Therefore, 9w > On even if any jet break in the optical 
occurs slightly after the jet break in the X-ray. This case 
can be applied, for example, to GRB060605. Thus, we 
conclude that Scenario A" fits better the cases of GRBs 
that show optical breaks only slightly later than the break 
in the X-ray. 

Note, though, that scenario A" can be solved even with high 
values of the kinetic energies. The following choice of param- 
eters satisfy the conditions; £'52,0,^ = 3 x 10^, ee,-i,n = 0.1, 

es,-2,n = 2 X 10 £52,0,11) = 200, ee,_l,iu = 0.1, £3,-2, id, 

n = 0.5. 

A possible advantage of this scenario is that it does not 
necessarily require high values of kinetic energy of the 
ejecta, so it can be applied to dim bursts and/or bursts 
with higher efficiency rj with respect to other models. 

As a potential drawback, in Scenario A" fine tuning is 
not removed, because a few inequalities are satisfied within 
factors of 1.5-2. 
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fx.w < -.fx,n 



4.4 Scenario B 

In this scenario, the conditions to be fulfilled are: 
fo,^ > 2/o,n (61) 

(62) 

i^x > T^c,w at 300 s after the trigger (63) 
Vm.m < i^o at 300 s after the trigger. (64) 
Vc,m > at 0.1 d after the trigger, (65) 

i^x > Vc,n at 300 s after the trigger (66) 
I'x > J^m,n at 300 s after the trigger (67) 
i^m,n > i^o at 0. Id after the trigger (68) 
which now give 

(69) 



0.85 pi. 35 1.4 „Q t:,l/3 p5/6 -3 



0.1 j^l.l 1.4 , 1 0.1 r-,1.1 1.4 

^B,-2,w-t^52,0 ,n^e, — 1 ,ri 7 

-3/2 p-1/2 -1 ^ 09 c 

<=S,-2,ti.-'^52,0,ii.'^ <^^.0, 



(70) 
(71) 
(72) 
(73) 
(74) 
(75) 
(76) 



1/2 2 pl/2 „ „ 

'=S,-2.m'=e,-l,™-'^52,0,tu ^ ^■'^ ' 

^fl'^2,»^^5"2ll'^"' >0.22. 

^fl'^2,„^5"2!o^>"' <32.5, 

<^S,-2,ueE,-l,n-E5^,0,n < 4 X 10 , 

1/2 2 ril/2 ^ 1 /I 7 

<=S,-2,n<=e,-l,n-e'52,0,7i > 1*. |- , 

By comparing Eq. [69] and 1701 we can immediately in- 
fer that ee,-i,n should be as high as possible, in order for 
these two equations to be fulfilled more easily. Furthermore, 
a high value of ee,-i,n makes easier to have the synchrotron 
frequency of the narrow component higher than the opti- 
cal band (Eq. I76[l even at late times. Therefore, we assume 
again the equipartition value of ee,-i,n = 3.3, as in the pre- 
vious scenarios. 

It is easy to check that Eq. [76] gives limits less stringent 
than Eqs. I74I75I By combining the latter two, we can see 
that they are satisfied by all plausible values of the energy 
of the narrow component (since they only imply i?52,o,n < 
x7.3 X lO^n), while for the other parameters they give 

e'"-'^ < 3£;'/'' n'/^ 

'^B,-2,n ^ ''-^52,0,71 i 

eR'_o „n-' < 1.2 X 10' . 



(77) 
(78) 



Therefore, as far as the narrow component is concerned, 
scenario B requires that n and eB,-2,n are not simultane- 
ously very small. For instance, for i?52,o,n ~ 10, it must be 
es,-2,n > 5 X 10~^n~^''' and for value of values of es,-2,n 
the corresponding limit on n must be computed accounting 
for Eq. [7S] as well. 

Stringent limits on the eB,-2,w can be obtained by con- 
sidering the conditions of the wide component. By using 
Eq.EIland El we have 



-^52,0,11. < 0.35neg__i_, 



(79) 



As we can see, ee,_i^^ should be quite small, in order to 
permit values of kinetic energy of the wide outfiow that are 



comparable with those observed in a few luminous GRBs, 
of the order « 10^'* ergs (Frail et al. 2001). For instance, 
if n = 0.1 and £52,0,™ = 250, then ee,-\,w < 0.2. In the 
following we will assume ee,--i,w = 0.06. Finally, an upper 
limit on £3,-2,111 can be then obtained from Eg. 1721 es,_2,uj < 

0. 05eJ'ij ni-^52 w which, with our choice of the parameter 
values, requires es,_2,™ < 40. 

It can be easily shown that, by using Es2,o,n = 30, 
n — 0.005, es,-2,n = 1-5, ee,-i,n = 3.3, £52,0,^ = 300, 
es,-2,iD ~ 10, ee,-i,w ~ 0.06 all the required conditions are 
satisfied. We notice that this scenario again requires a large 
degree of fine tuning between the parameters. It also requires 
a high value of kinetic energy of the wide component, almost 
as high as in Scenario A'. It therefore requires that the ef- 
ficiency of conversion of this kinetic energy into 7-rays is as 
low as in A'. Furthermore, in Scenario B, the relative ratio of 
the two component isotropic energies is Ez2,n/ Ez2,w ~ 10%, 

1. e. lower than in Scenario A'. Such a large difference in 
the two energies might cause the beaming angles of the two 
outflow not to differ considerably, unless a jet break in the 
optical occurs much later than the break in the X-ray (see 
Eq.IIJ. 

4.5 Scenario B' 

We will now explore a variant of the previous case, in which 
the cooling frequency of the wide component lie above the X- 
ray band. We therefore reverse condition 1631 Notice that the 
time when this condition has to hold changes as well; it can 
be shown that, in this scenario, if it holds at 0.1 days then 
it also holds since the beginning. Note also that expression 
1701 that relates to condition 1621 has to be changed as well. 
Overall, the required conditions now read; 



fo,w > 2/0,71 






(80) 


1 J- 

JX.w < -JX,n 




(81) 




at 300 s after the trigg 


;er 


(82) 


J^x ^ ^c,w 


at 0.1 d after the trigf 


rer 


(83) 


^x ^ ^m.n 


at 300 s after the trig 


ger 


(84) 


Vm,n > VO 


at 0. 1 d after the trig 


;ger 


(85) 


l'm,w < fO 


at 300 s after the tri; 




(86) 



which translate into: 
1.35 1.4 

52,0,tii^e,-l.tu 



0.85 r,1.35 1.4 ^ on c; 1/3 

^s, -2. 111^^52. n. -1 .1/, > /y.oe 



p5/6 -2/3 
2,n. 52,0,n'^e,-l,n ' 



(87) 



0.85 T^l.35 1.4 1/2 
£s,-2,ii>.C'52,0 

-■A/2 -1/2 1 

^B,-2.n^52,0,n'' <kOZ.O, 



^ o c^^in-2 0.1 ril.l 1.4 
< Z.OXiU eB,-2,n-C'52,0 



-3/2 „-l/2 „-l ^ 4 X in2 
'=S,-2,ii,-'^52,0,tii"' > 4 X iU 



^B'-2.n<-l,nEt^2:0,n < 4 X 10^ 



1/2 v?^/ ^ \ 1 /I 7 

'=S,-2,n<=e,-l,n-e'52,0,7i > 1*. ^ , 

1/2 



^1/2 
^52,0, 

^1/2 
52,0, 

pl/2 



■=S,-2,m"=e,-l,™-E52,0,tu < 0-22 . 



(89) 
(90) 
(91) 
(92) 
(93) 



As it can be easily seen, the simultaneous validity of 
both ineaualities l87l and [88](which have similar left members 
apart from the factor n^^^) crucially depends on the value 
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of £e,_i,„, that must be relatively large. Therefore in the 
following we assume again ee,_i,„ = 3.3. Once ee,-i,n has 
been assigned, Eq. |87] and |88] are more easily satisfied for 
relatively low values of the density and of eB,-2,n and for 
relatively high values of the kinetic energy of the narrow 
component. Besides, since ee,-i,n = 3.3, relations [77] and 
1781 involving the narrow component only, still apply. Based 
on that, we assume the following set of parameters for the 
narrow component: i?52,o,n = 4000, ee,-i,n = 3.3, es,-2,n = 
0.2 and a density n = lO"'^. With these choices, some of the 
Eqs. I87l l 92l are trivially satisfied, while the others give 

(&^,-2,wEl2,0,ni<^lf-\.w > 7.5 X 10^ , (94) 

eB^.^ElJl^el-t^.^n"/-' < 1.05 x 10^ , (95) 
^~B%^E;,]il>4. (96) 

From Eq. [S3] we can isolate an expression for ee,-i,w 
which, substituted into Eg 1941 gives 

£3^2,^,^52,0,™ > 2.3 X lO-* . (97) 

From this last equation we can immediately infer a lower 
limit on the value of £52,0,™ ■ Since the highest theoretical 
value of eB,-2,m is 33, achieved at equipartition, the mini- 
mum value of E^2,Q,w = 4 x 10^, which is admittedly very 
high. By using this value of i?52,o,u, in Eq. [96] we derive 
an upper limit on £5,-2,™ < 6 x 10"'^. Also, we can obtain 
an upper limit on ee,-i,™ by solving Eq. [^for eB,-2,w and 
substituting the resulting expression into Eq. 1931 We obtain: 

ee,-i,»£;5"2!o,l < 9.6 X 10-^ . (98) 

By using the upper limit on £52,0.11, quoted above, this 
last equation gives ee,-i,ii, < 2.5 x lO"'^. It is easy to verify 
that, for these values of the parameters of the wide outflow, 
Eq. [M] cannot be satisfled, unless £52,0,11, is unphysically 
large, ~ 10^. Scenario B' therefore cannot be assumed in 
our model. 

4.6 Scenario B" 

We will now explore a variant of scenario B, in which the 
synchrotron peak frequencies of both components are below 
the optical band, and the cooling frequencies are between the 
optical and the X-ray band. Overall, the required conditions 
now read: 



/o,ii, > 2/0, n 






(99) 


JX,w < -JX,n 






(100) 


Vm.w < VO 


at 300 s after the trif 


;ger, 


(101) 


^ ^C,lt, 


at 300 s after the trigj 


;er 


(102) 


VO < Vc,n} 


at 0. 1 d after the trig; 


jer 


(103) 


VO > I^m,n 


at 300 s after the trig 


;ger 


(104) 


Vc,n < l^X 


at 300 s after the trig] 


ger 


(105) 


Vc,n > VO 


at 0.1 d after the trigf 


;er. 


(106) 



which translate into: 

0.85 i;,1.35 1.4 . „ 0.85 r7,l-35 1.4 /in'7\ 
es,-2,ii,-C'52,0 > ^efl,-2,n-C'52,0 ,n^e, — l,ri i i ) 



0.1 z^-'---^ ^ 
£s,~2,ii,-C'52,0,ii,£e, 


,4 1 0.1 1-4 
,-1,111 > 2^B,-2,n-t^52,0,n^e,-l,n ) 


(108) 


1/2 ^1/2 2 
'=3,-2,11,-^52,0,11,'^=: 


,-i,i„ < 0.22 


(109) 


-3/2^-1/2 -1 


> 0.22 


(110) 


-3/2 „-l/2 - 
'=3,-2,11,-^52,0,11," 


< 32.5 


(111) 


1/2 „l/2 2 
'=S,-2,n-^52,0,n'=e, 


-l,n < 0.22 


(112) 


-3/2 ^,-1/2 - 
'=S, -2,11-^52,0, 


> 0.22 


(113) 


-3/2 ci-1/2 - 
^S, -2, 11-^^52, 0,n'^ 


^ < 32.5 


(114) 



Placing a different condition fo,m > 2fo,n and because 
v,n,n is below the optical band since the very beginning, it 
is no longer necessary to assume a high value for ee,-i,n to 
simultaneously satisfy the flrst two conditions. Instead, by 
combining Ea lll2l and lll3l we obtain 

£52,0,11 < 3.5e;r^i^„n (115) 

From this inequality, we derive that ee,_i,„ should be 
small, to allow high values of kinetic energy of the narrow 
outflow. As for the wide outflow, condition Eg 1791 still 
applies. 

For the following values of parameters, all the relevant 
inequalities of scenario B" are satisfied: £52,0,11, = 0.25, 
ee,-i,ii, = 0.25, es,-2,ii, = 10, £52,0,11 = 0.5, ee,-i,n = 0.25, 
es,-2,n ~ 0.3, n — 0.75. Scenario B" can be solved for 
higher values of kinetic energies as well: £52,0,11, = 20, 
^e,-i,w ~ 0.1, es,-2,ii, ~ 5, £52,0,11 ~ 90, ee,_i,„ = 0.1, 
es,-2,n ~ 0.15, n = 0.75 satisfy all conditions. 

Scenario B" is similar to Scenario A" , in the sense that 
it can be resolved for high and low values of the kinetic 
energies, and even in this case, £52,0,11, < £52,0,11. Likewise, 
Scenario B" does not solve the problem of fine tuning, and 
the ratio of £52,0,11, /£52,o,n is is much higher than in A". 
This scenario cannot thus be employed for cases in which a 
jet break occurs in the optical slightly after the jet break in 
the X-ray. Furthermore, it still presents the problem of fine 
tuning. 

4.7 Summary 

In summary, we have shown that there are at least two sce- 
narios of "A" kind and two of "B" kind that are satisfied for 
non-unreasonable values of the parameters. A drawback is 
that in all cases we require a large degree of fine tuning, since 
the allowed region in the parameter space is small. Since the 
bursts with chromatic breaks may not be rare (Liang et al. 
2008), fine tuning can represent a problem for our model. 

We would like now to address the point of the reliabil- 
ity of our model at late times, i.e. after the break observed 
at 0.1 days after the trigger. Within our model, this break 
is interpreted as a jet-break. This implies that, from this 
time onwards, the fiux of the narrow component is expected 
to decrease considerably faster than before, while the flux 
due to the wide outflow does not change its decay slope. 
Therefore, it is important to check that the flux in the X- 
ray due to the narrow component remains above that due 
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to the wide component even at late times. Would this con- 
dition not be satisfied we should observe a flattening of the 
X-ray lightcurve, as fx,w becomes comparable to fx,n at 
some time after the end of the plateau phase; this is clearly 
not observed in our GRB liglitcurvcs. 

Now, for p = 2.4, the X-ray flux of the narrow com- 
ponent decreases with time as fx,n oc t~^'^ in Scenario A' 
and A", and fx,n oc t~^'^'' in B and B" respectively. fx,w 
always decays as ~ t~° ''^. Therefore, the ratio fx,n/fx,w 
will decrease as t^^''^^ in scenario A' and A" and as in 
scenario B and B". With our suggested choice of parameters, 
condition fx,n > fx,w is satisfied (by a factor of ~ 10) in 
scenario A' at 0.1 d after the trigger, suggesting that a flat- 
tening of the X-ray lightcurve will not be seen before 2-2.5 
days after the trigger, when lightcurves are usually poorly 
sampled. In scenario A", fx,n > fx,w by a factor ~ 2.5 only 
at 0.1 d after the trigger; therefore, in this case, the X-ray 
flux of the wide component becomes comparable with that 
produced by the narrow outflow as early as ~ 0.3 d after the 
trigger and the X-ray decay slope should become similar to 
that in the optical, unless an early jet break occurs in the 
wide outflow as well. 

In the case of Scenario B, fx,n > .fx,w is satisflcd by a 
factor of ~ 20 0.1 d after the trigger. One should thus expect 
a flattening as late as in Scenario A". Finally, in Scenario 
B", fx,n ~ 2.5/x,n; for our choice of parameters, therefore 
the same restrictions of Scenario A" apply in this case, too. 

In drawing our scenario, we restricted ourselves to the 
simplest case of side-spreading jets and a constant density 
medium, with the addition of energy injection (Panaitescu 
ct al. 2006b) parameterized as L oc t'^'^ , and flxed p = 2.4. 
We also assumed a simple hierarchy between the relevant 
frequencies. In this simplified case, we have shown that our 
model successfully explains the characteristics of all bursts 
in our sample, with the only difl^erence that in some CEises we 
need to assume vx > Vc, and in some others the reversed in- 
equality. In many cases, the fraction of energy of the narrow 
outflow given to the emitting electrons has to be close to the 
maximum value allowed for adiabatic expansion (Freedman 
& Waxman 2001, Yost et al. 2003). In the case of GRBs 
without well sampled optical emission, we have deemed not 
to assume the Scenarios A" and B", which would require 
the presence of a flattening of the X-ray lightcurve only a 
fraction of day after the trigger, which is not observed. 

It is worth mentioning that we have also explored Sce- 
narios A', B, A" and B" in a wind scenario. We adopted 
the same frequency hierarchies of these two cases, but we 
replaced the set of equations with the set that describes the 
evolution of the characteristics frequencies and peak flux 
in a wind environment. These formulae were taken from 
Yost ct al. 2003. We found that, even in the case of cir- 
cumburst medium environment, these four scenarios basi- 
cally reproduce the observed behaviours, but flne tuning is 
not removed. Of course, it is possible to apply more compli- 
cated scenarios. For example, we may choose values of the 
parameters q and p which are different from those we have 
adopted in this paper, to reflect intrinsic differences among 
the various bursts. Changes of p and q from the values we 
have taken would result in a modification of both the expo- 
nents and the coefficients of the mathematical expressions 
we have used so far. As a consequence, some scenarios might 
not be viable anymore, or others could become applicable. 



We notice that our model can easily explain one of the 
most striking characteristics of the GRBs studied by Swift, 
i.e. the lack of evident jet breaks in the X-ray lightcurves 
(Burrows & Racusin 2007). In our scenario jet breaks are 
actually observed, but they arc not so steep as wo would 
expect from the traditional closure relationships (Sari, Pi- 
ran & Helpern 1999) due to the ongoing energy injection. 
Our model predicts that the steep decay slopes, like those 
observed in the optical in pre-Swift GRBs at late times, 
arc possible only once the energy injection has terminated. 
We note that our model might, in principle, be extended to 
all GRBs featuring the canonical lightcurve (Nousek et al. 
2006), even those without chromatic breaks. The implica- 
tion would be that, in those cases where optical and X-rays 
lightcurves show a simultaneous break at the end of the slow 
decline phase, the emission in both bands would arise from 
the same outflow. However, in our scenario the break is not 
caused by the end of an energy injection phase, as generally 
assumed when interpreting the canonical lightcurve, but by 
a jet break. Once the energy injection has terminated, the 
decline slope of optical and X-ray fluxes will assume the 
more typical values of a ~ 2. Thus our model can also ex- 
plain GRBs which show achromatic breaks only. The values 
of the decay and spectral slopes of the GRBs we have studied 
in this paper are not uncommon, supporting the idea that 
our model could be applied in several cases. Our interpre- 
tation can call for a deep revision of GRB physics, such as 
the mechanism that produces the outflow and the energetics 
involved in the process. We need to explain how the cen- 
tral engine can either be active for several days, or produce 
a long trail of shells that merge for such a long time. Be- 
sides, we should flnd mechanisms that can commonly beam 
ejccta into cones, which can have opening angles as narrow 
as 5 X 10~^ rad. 



5 CONCLUSIONS 

In this paper, we have reanalysed the full sample of 
Swift GRBs whith chromatic breaks, originally discussed by 
Panaitescu ct al. (2006a). In addition, we have also stud- 
ied GRB 060605, another Swift burst with good quality 
XRT and UVOT data and a chromatic break in the XRT 
lightcurve. We have shown how our model, based on a pro- 
longed energy injection into a double component outflow 
and a jet break, is physically plausible and can well ex- 
plain the behaviour of the optical and X-ray emission of 
GRB050319, 060505 and GRB050802 (see also Gates et al. 
2007). GRB050922c has been shown not to require a chro- 
matic break. We note that our model can also be applied 
to the other GRBs with claim of chromatic breaks pub- 
lished in Panaitescu et al. (2006a) and might, in principle, 
be extended to all GRBs featuring the canonical lightcurve 
(Nousek ct al. 2006), even those without chromatic breaks. 
We emphasize that it would have not been possible to de- 
rive our conclusions if we had considered the X-ray data 
only, since GRB050319, GRB060605 and GRB050802 ex- 
hibit a canonical X-ray lightcurve. Instead, the combined 
optical aiid X-ray analysis has shown that the component 
responsible for the optical is uncoupled from the outflow 
that produces the X-ray emission. In our model, the ejecta 
responsible for the X-ray emission are narrowly beamed, and 
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undergo an early jet break that explains the chromatic break 
seen in the X-ray only. Our model of combined jet expan- 
sion and energy injection may have deep consequences on 
our understanding of the GRB, since it calls for a revision 
of the physics processes that take place in these objects. 
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Figure 1. Lightcurvcs and SEDs of GRB050319, GRB050802, GRB060605. UVOT lightcurves are fitted by simple powerlaws, while 
XRT lightcurves and SEDs are fitted by broken powerlaws 
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Figure 3. The 6 configurations of the wide (W) and narrow (N) components explored in the text. Vertical lines indicate the optical and 
X-ray band respectively. 



